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In an attempt to understand quantitatively the remarkable 
discoveries of metal-insulator transitions in two-dimensional 
systems, we generalize Mott's variable range hopping theory 
to the situation with strong Coulomb interaction. In our for- 
mulation, the Gaussian form is adopted into the expression of 
the hopping probability, and the effect of Coulomb gap is also 
considered. After taking account of the newly proposed scal- 
ing consideration, we produce the dynamical and localization 
length exponents, which are consistent with the experiments. 
We then clarify the physical content of our formulation and 
explain the universality of the localization length exponent 
suggested by a series of experiments. We also discuss the 
general scaling function of both temperature and electrical 
field on the insulating side of the transition. 

The conventional scaling theory for disordered systems 
holds that without interactions any extent of disorder is 
sufficient to localize the electrons and no true metallic 
behavior is possible at T = in two dimensions JjJ . The 
inclusion of interactions to this picture proved to be much 
more difficult, and no satisfactory theoretical picture has 
emerged in spite of decades of continuous effort. However 
the belief that all the states are localized at d = 2 has 
remained unchallenged. 

Recently, a series of surprising experiments have been 
conducted on 2D electron gas and hole gas in zero mag- 
netic field, which demonstrated convincing evidences in 
support of a true metal-insulator transition(MIT) in 2D 
|| . In these experiments, The resistivity scales with tem- 
perature or electrical field with a single scaling parameter 
that approaches zero at a critical carrier density n c . Fur- 
ther more the scaling function of resistivity near the tran- 
sition possesses a very simple exponential form, which 
leads to the reflection symmetry across the transition, as 
emphasized in ref H- That is, 



p(Sn,T) = l/p(-6n,T), 



(1) 



where 5n = (n — n c )/n c , p is the resistivity measured 
in unites of its critical value p c . This symmetry is very 
similar to some recent experimental results in quantum 
Hall systems although the mechanism may be quite 
different || . These remarkable discoveries have evoked a 
wave of theoretical efforts to understand them || M . Al- 
though it is generally agreed that strong interactions play 
an important role, no specific microscopic mechanism has 
been able to explain the experiments both qualitatively 
and quantitatively. 



In a recent Letter pj], Dobrosavljevic et al. gave a 
general scaling consideration of the 2D MITs. In their 
formulation, they treated the temperature-dependence of 
the conductance in the quantum critical region. After 
assuming the MIT occurs at g = g c where f3{g c ) = 
and introducing the scaling variable t = \og(g/g c ), they 
considered the linear approximation to /3(f) near t = 
and obtained 



/3(f) 



dt 



t 



0(f 2 ). 



d(logL) v 
By integrating from I to L, they got 

where to = log(go/g c ) is given by the value go of the 
conductance at the microscopic scale I. Then after using 
fo ~ (go — 9c) I gc oc Sn , where Sn = (n — n c )/n c and n 
is the density, they reached 

g{L) = g c exp[A5n(L/l)$]. 

Then representing the length scale L by the temperature 
as T oc L~ z , one finally arrived at 



g(8n,T) = g c e~xp(A6n/T "* ). 



(2) 



This expression gives the explicit exponential form of the 
scaling function of the conductance with respect to tem- 
perature except for the exponents v and z to be deter- 
mined, and the reflection symmetry in eq([j]) follows from 
it directly, which are precisely what has been seen in 
the experiments. As emphasized by the authors JtJ, the 
above relation is only valid inside the critical region, or 
for T > Tq on \5n\ l/z , which is also verified experimen- 
tally S. In what follows we will use this result as a re- 
strictive condition to our general derivation, which is the 
generalization of Mott's well-known variable range hop- 
ping (VRH) theory. Then we will arrive at the explicit 
scaling function of the resistivity, from which the expo- 
nents v = 1.5 and z = 1 can be read off easily, which are 
consistent with the experiments || . We then discuss the 
physical mechanism of our formulation and explain the 
universality of the exponent v = 1.5 Q, which is charac- 
teristic of a general two dimensional systems with strong, 
long-range interactions( like Coulomb interaction). At 
last, we discuss the general scaling function of both tem- 
perature and electrical held. 
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Before starting with our formulation, let us give an 
introduction to some background that motivates the 
present work. In an influential work ]ic[ ], Mott suggested 
his celebrated theory of VRH which described an impor- 
tant conducting mechanism in a localized electron sys- 
tem. His central idea is to optimize the hopping proba- 
bility with respect to the hopping distance R, which gave 
the temperature-dependence of the conductivity 

a(T) oc exp[-(^)^r], 

where the dimension d = 2 here. This result has received 
substantial experimental supports in the past decades 
(□J. Later, Efros et al. made an attempt to include 
the effect of Coulomb interaction into the VRH theory 
They only considered the Coulomb gap correction 
to the density of states and obtained 

CT (T) oc exp[-(^)5], 

This has also been verified in the experiments performed 
on dilute electron systems where Coulomb interaction 
dominates the transport properties |fj"3f . Although this 
relation can produce the right exponent z — 1, it is never- 
theless inconsistent with the concrete expressions derived 
from the accurately designed measurements performed 
recently |2j and fails to explain the general localization 
length exponent v w 1.5 revealed by more than one ex- 
periment. 

In the following, we will extend the VRH idea to the 
situation with strong Coulomb interaction. Our general- 
ization goes beyond the one attempted by Efros et al., 
and aims at a better understanding of the recently surg- 
ing discoveries of MITs in two dimensional systems, es- 
pecially the exponent v, which seems to be a universal 
exponent. 

The hopping probability for a particle to cover a dis- 
tance of R is 



p(R) oc exp(— a-^j 



^ 2 



(3) 



where £; is the localization length on the insulator side of 
the MIT, which diverges as £ oc | "^" c \~ w when approach- 
ing the MIT (y is the localization length exponent); a is 
a non-singular parameter, A is the energy difference be- 
tween the initial state and the final state for the hopping 
particle. The only difference from Mott's equation is the 
Gaussian form dependence of p{R) on R, which is vital 
to the final conclusion and contains a profound physics 
content (see below). If there is no interaction, we simply 
have 

A cx [R d N{E F )]-\ 

as it is supposed in Mott's derivation, where N(Ep) is 
the density of states at the Fermi level Ep. However, for 



the case of strong Coulomb interaction, the Coulomb gap 
will depress the density of states near the Fermi surface. 
Following Efros et al. |Q, we can write 

N(E) oc \E- E F \, 

which means N(Ep) — 0. Therefore N(Ep) in the above 
equation should be substituted by the following integra- 
tion which counts the number of states available for a 
R-range- hopping: 



N(E f 



AE 

oc AE 



1 f AE 

— / dEN(E + Ep) 
^ E Jo 



(4) 
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1 

R' 



where a linear dispersion AE oc Ak has been used, and 
Afc is the momentum quanta proportional to 1/R. 

Therefore A oc l/R d ~ l . This means the effect of 
Coulomb interaction is to reduce the dimension by 1. 
Therefore in what follows, we will use the effective di- 
mension d' = d — 1. 

Then we maximize cq(|J) with respect to R, and obtain 



£ 2 i 
oc (— )3H3. 



Therefore we finally arrive at 

(|An|)^+3 

cr(n,T) (xp(R max ) oc exp{ 5 }. 



(5) 



This is the central equation of the present paper, which 
should be contrasted with the relations of Mott and Efros 
et al.. 

Taking into account the classical scaling relation 
a(n,T) = F(^), 

where F(x) is the scaling function, b = 1/zv, and z is 
the dynamical exponent, we can easily see 

z = d' = 1 

for a two dimensional system. This result is well-known 
for a Coulomb interaction dominated system ]l5[| . This 
relation is also verified in many recent experiments on 
2D MIT @. 

However, there is more to obtain from eq(|5|). Since the 
above single parameter scaling argument has given the 
following restrictive form 

p{n,T) oc exp(-A^), 

we have an additional condition to satisfy: 
2d'v _ 1 



2 



which gives 



d! 



2d' 



1.5. 



This is very close to the localization length exponent re- 
ported widely by a series of experiments Therefore, 



JAnL 

p(n,T) oc exp(^^). 



(6) 



We would like to comment that even if we have adopted 
the traditional p oc exp(— R/£) instead of the Gaussian 
form, we can still arrive at z = 1. However the desirable 
v can not be obtained in that way. That is to say, the 
traditional VRH theory is not consistent with the scaling 
argument presented in rcf fi, under the condition that 
v = 1.5, z = 1. 

Now let us give some comments on the physics con- 
tained in the Gaussian form adopted in eq(|3j) . This form 
is quite natural and reasonable in physics, because of its 
direct connection with the eigenstates of the harmonic 
oscillator(HO), one of the few universal and analytically 
transparent models. Motions in a general smooth enough 
potential near its minimum can be approximately de- 
scribed by a HO. Therefore, if there does exist such ef- 
fective potential that can equivalently describe the mo- 
tion of the localized particles relevant to the 2D MITs, 
the Gaussian form will be probably a most suitable one 
to be employed. However, considering the success of 
Mott's original relation and the strong Coulomb inter- 
action present in the samples used for the research of 
2D MITs, we are more inclined to the viewpoint that 
the effective potential felt by the hopping particles are 
probably relevant to the long-range Coulomb interaction. 
First, the property of long range ensures that the poten- 
tial experienced by one particle is the sum of the contri- 
butions of many other particles, which is averaged to be 
smooth and slowly changing. Secondly, the domination 
of Coulomb interaction over the kinetic energy suppresses 
the dynamical fluctuations which push the particles away 
from their most stable position with the lowest potential 
energy. Therefore it will be a relatively good approxi- 
mation to represent the effective potential for a hopping 
particle by an HO potential, with the particle's initial 
position at the center. For comparison, we give further 
comments on the traditional form exp(— R/t;). This form 
is specific to a square well potential, which is widely used 
in models of localization. As we believe, the above form is 
good for short-range interactions like a hard-core poten- 
tial, or in a tight bound model where electrons are bound 
by the ionic potential. However, in case of long-range in- 
teractions, the much smoother HO potential is better. 
This argument may explain why v = 1.5 is extensively 
reported in strongly interacting 2D systems, independent 
of the sample parameters. [|l(| 

Then we discuss what happens when an electrical field 
E is turned on, and what will be the explicit form of the 



scaling function with two variables E and T on the insu- 
lating side of the MIT. With the assistance of electrical 
field E, the hopping probability can be written as 



I? 2 



A 



ER 



p(R) oc exp(-a- - ^ , ^ 



(7) 



With z = 1 in mind, we can define two dimensionlcss 
scaling variables: 



Ee 



t = T£. 



Then suppose R ~ £e v t x (x, y are exponents to be de- 
cided), and substitute it into eq(0), we have 

p(R) tx exp[F(e,t)], 

where F(e, t) is the general scaling function of e and t, 
and is the sum of the following three parts |0| : 



Fi ~ -e 2y t 2x , 



(8) 



F 3 ~ e»^ L t; L 

We then maximize the eq(Q) with respect to R, and get 



2aR 



i 



E 



R 2 T T ' 



By changing into e and t the above equation can be writ- 
ten as 

eV-e^r^'+ef 1 . (9) 

Therefore the general scaling function is the sum of Fi,F 2 
and F 3 , with the exponents x and y decided by the con- 
striction equation (||). This suggests that the scaling 
function will be of different forms for different regions 
of e and t, and demonstrates complicated behaviors as a 
result. 

We first discuss the simple case of |e| << 1. In order 
for F(e,t) to give a convergent result, we must satisfy: 

2y>0, -y>0, y + l>0. 

Therefore it is necessary that y — 0, which is substituted 
into eq(||). Then we get x = —2x — 1, so x = —1/3. We 
thus arrive at the scaling function for |e| << 1, 



F(e,t) 

For e = 0, we obtain 



-t 



-2/3 



+ et 



-4/3 



p(n, T) oc exp[— F(0, T)} oc exp(- 



~ASn 

T2/3 



), 



which is exactly eq(|^) we get earlier. 

It is also interesting to understand the behavior of the 
scaling function for t << 1 and t >> 1 which will be 
reported elsewhere because of its complexity. 
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Before closing, we would like to give more comments 
on our use of a hopping theory to describe transport 
properties near MITs. We first point out an important 
difference of our hopping picture from Mott's original 
theory. In Mott's formulation, the relevant system is 
strongly localized, so the localization length £ is very 
small and the overlap between adjacent wave-packets is 
also small. Therefore, in order for a hopping precess to 
contribute to the conductance, the hopping distance R 
must be larger than £, which puts a constraint on R m ax- 
This leads to the conclusion that Mott's hopping is only 
valid at low enough temperature, or T < To. However 
the situation near MITs is quite different, where £ be- 
comes very large, and the overlap between adjacent wave 
packets is also large, so a particle hopping between neigh- 
boring sites can contribute to the conductance without 
covering a distance as long as £. Therefore the low- 
temperature-constraint on hopping mechanism is actu- 
ally absent, and we are justified to apply it inside the 
critical region T > Tq. The relevance of hopping for 
both T > T and T < T is supported by inspecting the 
inset of Fig. 2 in ref H, where no crossover behavior is de- 
tected near Tq for the insulating phase, while the metallic 
phase shows obvious deviation from the linear behavior 
for T <Tq. In fact, the use of hoppings to study critical 
phenomena is not new jTsI, and its relevance to critical 
exponents (i.e. z) was already established by Efros et al. 

In conclusion, we have presented a generalization of 
Mott's variable range hopping theory to the situation 
with strong Coulomb interaction in order to under- 
stand quantitatively the remarkable discoveries of metal- 
insulator transitions in two-dimensional systems. In this 
formulation, the Gaussian form is adopted into the ex- 
pression of the hopping probability, and the effect of 
Coulomb gap is also considered. After taking into ac- 
count the newly proposed scaling consideration, we for 
the first time give the explicit scaling function of the re- 
sistivity on the insulating side, and as a result determine 
the dynamical and localization length exponents, which 
are consistent with the experiments. We then clarify the 
physical meaning contained in our formulation with its 
relevance to the long-range Coulomb interaction. We also 
discuss the general scaling function of temperature and 
electrical field on the insulating side of the transition. 

We acknowledge helpful discussions with Z. B. Su and 
V. Dobrosavljevic. This work is in part supported by 
National Natural Science Foundation of China. 
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